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Multiscaling in superfluid turbulence: A shell-model study 
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We examine the multiscaling behavior of the normal- and superfluid-velocity structure functions 
in three-dimensional superfluid turbulence by using a shell model for the three-dimensional (3D) 
Hall-Vinen-Bekharevich-Khalatnikov (HVBK) equations. Our 3D-HVBK shell model is based on 
the Gledzer-Okhitani-Yamada (GOY) shell model. We examine the dependence of the multiscaling 
exponents on the normal-fluid fraction and the mutual-friction coefficients. Our extensive study of 
the 3D-HVBK shell model shows that the multiscaling behavior of the velocity structure functions 
in superfluid turbulence is more complicated than it is in fluid turbulence. 
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I. INTRODUCTION 

The characterization of energy spectra and velocity 
structure functions fl} occupies a central place in the 
elucidation of the statistical properties of turbulence, be 
it in fluids M, conducting fluids Bum, or superflu¬ 
ids mm- For example, in fluid turbulence, we of¬ 
ten use the longitudinal velocity v structure function 
S p (r) = ([5u(r)] p ), where Sv(r) = [v(x + r) — v(x)] • [r/r], 
which scales as S p (r) ~ A p , for r =| r | in the inertial 
range < r < L; viscous dissipation is significant be¬ 
low the dissipation length scale rjd\ and L is the large 
length scale at which energy is injected into the fluid. 
The exponents (' p , which characterize multiscaling, are 
nonlinear, monotone increasing functions of p 0; simple 
scaling is obtained if £ p depends linearly on p , as in the 
K41 phenomenological approach of Kolmogorov [3® 
that yields (, p 41 = p/3. 

Direct numerical simulations (DNSs) play an impor¬ 
tant role in studies of structure-function multiscaling in 
fluid turbulence [lj, 0, 0; such DNSs have achieved im¬ 
pressive spatial resolutions (see, e.g., Refs. 0 , 0 ). By 
contrast, DNS studies of superfluid turbulence, whether 
at the level of the Gross-Pitaevskii (GP) equation 
( (2ll I® and references therein) or via the Hall-Vinen- 
Bekharevich-Khalatnikov (HVBK) equations ([3, ®] 
and references therein), have only achieved modest spa¬ 
tial resolutions. Furthermore, the large number of pa¬ 
rameters in these equations, e.g., the mutual-friction co¬ 
efficients, the ratio of the normal-fluid density to the su¬ 
perfluid density, and the Reynolds number, pose a signif¬ 
icant challenge for systematic studies of the multiscaling 
of normal-fluid- and superfluid-velocity structure func- 
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tions. It has been suggested, therefore, that shell models 
for the three-dimensional (3D) HVBK equations [3® 
be used first to study such multiscaling in detail. 

Ever since their introduction in the early work 
of Obukhov ®|i Desnyansky and Novikov ®], and 
Gledzer, and Ohkitani and Yamada ®, ®] (henceforth 
GOY), shell models have played valuable roles in eluci¬ 
dating the multiscalin g prop erties of structure functions 
of fluid turbulence [Tl, uHlI® . Over the years, such shell 
models have been used to study magnetohydrodynamic 
(MHD) turbulence Hall-MHD turbulence ®- 

®, fluid turbulence with polymer additives ®, fluid 
turbulence in two dimensions [50[, fluid turbulence in di¬ 
mensions in between two and three [Hlj ]. turbulence in 
binary-fluid mixtures Pvj and in rotating systems [53j |. 
and, as we have mentioned above, turbulence in superflu- 
ids [3® . Shell models have also been used to initiate 
studies of the dynamic multiscaling of time-dependent 
structure functions [3®. 

We build on the shell-model studies of Refs. ®-® 
to explore the dependence of the multiscaling exponents 
here on the parameters of the 3D HVBK model. It has 
been noted in Ref. [® that, given current computational 
resources, a systematic study of this parameter depen¬ 
dence lies beyond the scope of a well-resolved DNS of 
the 3D HVBK equations; however, such a study is possi¬ 
ble if we use shell models for these equations. Our study 
extends the work of Refs. [3 - 1® by obtaining a variety 
of results, which we summarize, before we present the 
details of our study. 

Our study of the 3D-HVBK shell model shows that the 
multiscaling behavior of the shell-model counterparts of 
velocity structure functions in superfluid turbulence is 
more rich than that reported in Ref. ®. Our results 
agree with those of Ref. ®] qualitatively insofar as we 
find that, in the limits when the normal-fluid fraction 
is either small (pure superfluid) or large, the equal-time 
multiscaling exponents are close to their classical-fluid- 
turbulence values. In addition, we And that there are 
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two regions, with intermediate values of the normal-fluid 
fraction, in which the multiscaling exponents are larger 
than those observed for the classical-fluid-turbulence or 
even Kolmogorov’s 1941 (K41) predictions [lj; between 
these two regions there is a region in which the multiscal- 
ing exponents are close to their K41 values. We have also 
investigated the dependence of the multiscaling expo¬ 
nents on the mutual-friction coefficient, with equal pro¬ 
portions of superfluid and normal-fluid components; here, 
our results show that, for small (weak-coupling limit) 
and large (strong-coupling limit) values of the mutual- 
friction coefficient, the multiscaling exponents tend to 
their classical-fluid-turbulence values, whereas, in an in¬ 
termediate range, there are deviations from the classical- 
fluid-turbulence behavior; in particular, the multiscaling 
exponents are larger than their classical-fluid-turbulence 
counterparts for high-order structure functions (order 

P > 3). 

The remainder of this paper is organized as follows. 
In Sec. HI] we describe the shell model and the numerical 
methods we use. Section m is devoted to our results. 
We end with conclusions in Sec. ED 


II. MODELS AND NUMERICAL SIMULATIONS 


The GOY shell model, for the 3D Navier-Stokes equa¬ 
tion mm for a fluid, comprises the following ordinary 
differential equations (ODEs): 


equations 03. j/17' is 

= — + p s a\7T + (2a) 

Du n o 

Pn^- = - — V P -PnVVT + ¥ n mf + V r y 2 U n , (2b) 

Ut p J 

with Du 1 /Dt = du l /dt + u* • Vu ! , the incompressibility 
condition V • u* = 0, and the superscript i G (n,s) de¬ 
notes the normal fluid (n) or the superfluid (s); p, a , and 
T are the pressure, specific entropy, and temperature, 
respectively; p n (p s ) is the normal-fluid (superfluid) den¬ 
sity; v n is the kinematic viscosity of the normal fluid. 
The mutual-friction terms, which model the interaction 
between the normal and superfluid components, can be 
written as F s mf = ~(p n /p)i m f and F^ = (p 8 /p)f m f in 
Eqs. (Eal) and (Eal) . respectively, where 


fmf 


B u) s 
2 |w s | 


B' 

X (u.'s X U ns ) + -^-UJ S X u ns , 


( 3 ) 


with u ns = (u n — u s ) the slip velocity, and B and B' the 
coefficients of mutual friction. In most of our studies we 
set B' = 0, so f m f = — -j|w s |u ns , which is the Gorter- 
Mellink form [Hj]. 

We use the following shell model for the 3D HVBK 
equations; it is based on the GOY shell model for a 
fluid H - 


d 

dt 


+ v n k 
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here we label shells by the positive integers m, in a log¬ 
arithmically discretized Fourier space, with scalar wave 
numbers k m = koX m , where ko = 2 -4 and A = 2. The * 
denotes complex conjugation, v is the kinematic viscos¬ 
ity, and u m (k m ) are the complex, scalar, shell velocities. 
The coefficients a = 1, b = —5, c = —(1 — 5) are cho¬ 
sen to conserve the shell-model analogs of energy and 
helieity in the limit of vanishing viscosity and the ab¬ 
sence of external forcing; the standard value of S is 1/2; 
N is the total number of shells; and f m is the external 
forcing, which is used to drive the system into a turbu¬ 
lent state that is statistically steady. The logarithmic 
discretization of Fourier space allows us to achieve very 
high Reynolds numbers, even with a moderate number 
of shells. In the GOY-shcll-model equations, direct in¬ 
teractions are limited to the nearest- and next-nearest- 
neighbor shells; in contrast, if we write the Navier-Stokes 
equation in Fourier space, every Fourier mode of the ve¬ 
locity is directly coupled to every other Fourier mode. 

The simplest form of the incompressible, 3D HVBK 



NL[<J+i^ + /s 


( 5 ) 


where 


NL[li m ] — l\cik rn U rn -\-\U rn -\-2 + 

T ck rn —2U rn — l'Um—2\ • 


( 6 ) 


Here, as in the GOY model, we have a logarithmically 
discretized Fourier space with shell-m wave numbers 
k = fcoA m , where kg = 2~ 4 and A = 2, and kinematic vis¬ 
cosities v n and v s for the normal fluid and the superfluid, 
respectively; of course, v s must vanish in a superfluid 
but, in practical numerical simulations, v n u s > 0 for 
numerical stability. The normal and superfluid dynami¬ 
cal variables are, respectively, the complex, scalar, shell 
velocities u^(k m ) and u^(fc TO ); and /)/ and /)/ are the 
external forcing terms. The coefficients a = l, b = — 1/2, 
c = —1/2 are chosen to conserve the shell-model analogs 
of energy and helieity in the limit of vanishing viscos¬ 
ity and the absence of external forcing. The shell-model 
analogs of the mutual-friction terms, which models the 
interaction between the normal and the superfluid com¬ 
ponents, are 


F s = 


p n BXl 


1/2 


2 p 




( 7 ) 
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and 


n 

F-=~^-^{u n m -u s m ). (8) 

2 P 

The shell-model superfluid and normal-fluid enstrophies 
are, respectively, 

n. = £^XJ a ( 9 ) 

m=l 

and 

N I 

(10) 

m—1 

The total energy is 

1 * 

E r =E n + E s = -^(| u «| 2 + |^| 2 ), (ii) 

m —1 

where -E„ and E s are the normal-fluid and superfluid 
energies, respectively. Other statistical quantities that 
we use in our study are as follows: The helicity is 


Hi = 



the energy spectra are 




1 I u. 


i 12 


2 k„ 


the root-mean-square velocities are 


(Ei 


|2 


1/2 


( 12 ) 


(13) 


(14) 


the Taylor microscale is 


The equal-time, order-p structure functions for the 
shell model are 


s;(k m ) = (K(t)<(f)] p/2 ) - kj ; , (i9) 


where the power-law dependence is obtained only if 
lies in the inertial range. The structure functions defined 
above show period-three oscillations because of three cy¬ 
cles in the static solutions of the GOY model for the 
Navier-Stokes equation 13311. Therefore, we use the mod¬ 
ified structure functions |3a l33j 


y l = 
^v ~ 


u m+2 u rn+l u rn ^ u m-l u m u m+l\ 


p/3 


r-C 


( 20 ) 

which filter out these oscillations effectively. The Sabra- 
model variant [25|, [26| of the 3D HVBK equations does 
not show such oscillations. We expect that the mul¬ 
tiscaling exponents Q, i £ ( n,s ), satisfy the follow¬ 
ing convexity inequality for any three positive integers 
Pi < P2 < P3 Q: 

(P3 ~ Pl)C l 2p 2 > (P3 - P2)C 2pi + (P2 ~ Pl)C 2 p 3 - (21) 


We obtain smooth energy spectra, without period-3 os¬ 
cillations, by using Ei(k m ) = T,\{k m )/k m , i £ ( n,s). 

To obtain a turbulent, but statistically steady, state, 
we force both the superfluid and the normal-fluid com¬ 
ponents with the forces 


/m’ s = (1 + *) x 5 x 10 -3 d ljm , (22) 


where <5i jm is the Kronecker delta. We use the second- 
order, slaved Adams-Bashforth scheme to integrate the 
3D-HVBK-shell-model Eqs. (J3J and © j32|,[59}- To study 
the multiscaling behaviors of structure functions here, we 
design the following three sets of runs: 

1. Gla-G9: In these runs, we use the values of p n /p 
and B , which have been measured at different tem¬ 
peratures in experiments on helium II [6Cj. We use 
suitable values of u n and v si which we list, along 
with other parameters, in Table HI 


\ _ Em F l (k m ) 

[Emk 2 m E*(k m )\ ’ 

the Taylor-microscale Reynolds number is 

FL&\ — Rrms A / y ■ 

the integral length scale is 

n _ Em F l {km)/km 

1 E ra^ikm) ’ 

and the large-eddy-turnover time is 

rpi _ 

eddy ~ h[x 4 5 


(15) 


(16) 


(17) 


(18) 


2. B1-B19: We vary p„/p between 0.05 —0.95 and keep 
B = 1.5 fixed. 

3. R1-R12: We vary B between 0.1 — 10 and keep 
Pn/P = 0.5 fixed. 

In the runs B1-B19 and R1-R12, we use v n = 10 -7 , 
v s = 10 -9 , and the time step At = 10~ 5 . 

We use the initial condition u^ s = (1 + i)k m e~ km , for 
1 < m < N, in the runs Gla-G9, PG1, and PG2; the GOY- 
shell-model runs PG1 (v n = 10~ 7 ) and PG2 {y n = 10 -9 ) 
are included for the purpose of comparison with the runs 
Gla-G9. In the runs B1-B19 and R1-R12, we use the initial 
values = u^ s k]l 2 e~ km e ldm , for 1 < m < N, where 
■dm is a random phase distributed uniformly on [0,27 t). 
We use the boundary conditions u l _ 2 = u l _ 1 = u l 0 = 0 
and u l N+1 = u l N+2 = 0, i £ (n,s). We report results 
for N = 36 shells; Ref. [24[ uses N = 18 and Ref. [26| 
presents data with N = 36. 


here and henceforth i £ (n, s). 
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Pn/p 

B 

Vn 

Vs 

At 

PG1 

- 

- 

io- 7 

- 

IO" 5 

PG2 

- 

— 

10 -9 

- 

10 -5 

Gla 

0.0450 

1.5260 

io- 7 

10“ 10 5.0 x 10 -6 

Gl 

0.0450 

1.5260 

io- 7 

10 -9 

10 -5 

G2 

0.0998 

1.3255 

io- 7 

10 -9 

10 -5 

G3 

0.2503 

1.0765 

10- 7 

IO" 9 

IO" 5 

G4 

0.4004 

0.9838 

io- 7 

10 -9 

10 -5 

G5 

0.4994 

0.9848 

io- 7 

10 -9 

10 -5 

G6 

0.6003 

1.0447 

io- 7 

10 -9 

10 -5 

G7 

0.6493 

1.1034 

io- 7 

10 -9 

10 -5 

G8 

0.6995 

1.1924 

io- 7 

10 -9 

10 -5 

G9 

0.7501 

1.3267 

io- 7 

10 -9 

10 -5 


TABLE I. Parameters for our 3D-shell-model runs (classical- 
fluid-turbulence) PG1, PG2 and 3D-HVBK-shcll-model runs 
Gla-G9: p n /p is the normal-fluid density fraction; B is the 
mutual-friction coefficient; v n ( v s ) is the normal-fluid (super- 
fluid) viscosity; At is the time step; we use N = 36 shells in 
our simulations. 


III. RESULTS 

We now present the results of our study of super¬ 
fluid and normal-fluid turbulence in the 3D-HVBK shell- 
model. We begin with energy spectra and then examine 
the parameter dependence of the exponents that charac¬ 
terize the multiscaling of structure functions. 

In Table HI! we list the values of A,, Re \, and T l eddy that 
we obtain from our 3D-HVBK-shell-model simulations 
PG1, PG2, and Gla-G9. Figure [Q compares E n (k m ) (full 
curves) and E s (k m ) (dashed curves) for four representa¬ 
tive values of p n /p (runs G1 (purple curves), G2 (green 
curves), G5 (sky-blue curves), and G9 (brown curves)). 
The inertial ranges of E n {k m ) and E s {k m ) exhibit scaling 
that is consistent with a fc -5 / 3 power-law form (orange, 
dashed line); of course, this exponent is not exactly —5/3 
if the structure functions display multiscaling. The runs 
PG1 and PG2 can be regarded as uncoupled (B = 0) nor¬ 
mal fluid and superfluid, respectively; we use them for 
the sake of comparison with other runs to show how the 
mutual friction modifies the energy spectra. When we 
couple the normal and superfluid components, as in the 
run Gl, E n {k m ) is pulled up towards E s (k m ), by virtue of 
the mutual-friction-induced tendency of locking between 
u n and u s (see Ref. E3) ; in contrast, in the absence of 
coupling, the spectra E(k m ) for the runs PG1 (yellow, 
full curves) and PG2 (yellow, dashed curves) lie far apart, 
especially in the dissipation range. 

We study the multiscaling behaviors of the velocity 
structure functions for the 3D-HVBK shell-model by cal¬ 
culating the multiscaling exponents and £®, for the 
normal fluid and superfluid components, respectively, by 
using the Eqs. (TSUI) for E 1 . In Table Hill in the Supplemen¬ 
tal Material, we list the values of these exponents, which 



A n 

A s 

Re"(xl0 6 ) 

Re 3 x (x 10 8 ) 

rjiTL 

eddy 

rjiS 

eddy 

PG1 

0.95 

- 

7.2 

- 

14.50 

- 

PG2 

0.50 

- 

310 

- 

20.33 

- 

Gla 

0.42 

0.28 

2.3 

15 

45.61 

45.61 

Gl 

0.70 

0.51 

4.2 

3.1 

21.80 

21.80 

G2 

0.73 

0.54 

4.4 

3.2 

22.03 

22.03 

G3 

0.82 

0.61 

5.2 

3.9 

19.83 

19.83 

G4 

0.71 

0.54 

4.4 

3.4 

18.41 

18.41 

G5 

0.89 

0.70 

6.2 

4.9 

17.15 

17.15 

G6 

0.94 

0.77 

6.9 

5.6 

15.45 

15.45 

G7 

0.94 

0.78 

7.0 

5.9 

14.94 

14.94 

G8 

0.95 

0.80 

7.2 

6.1 

14.51 

14.51 

G9 

0.95 

0.82 

7.3 

6.3 

14.42 

14.42 


TABLE II. Parameters from our shell-model runs PG1, PG2, 
and Gla-G9: X„ (A s ) is the Taylor microscale for the 
normal-fluid (superfluid); Rex (Re\) is the Taylor-microscale 
Reynolds number for the normal-fluid (superfluid); T addy 
{Teddy) i s the large-eddy-turnover time for the normal-fluid 
(superfluid). 



FIG. 1. Log-log (base 10) plots of the spectra E n (fc m ) (full 
curves) and E a {k m ) (dashed curves) from our shell-model 
runs: PG1 and PG2 (yellow curves); Gl (purple curves); G2 
(green curves); G5 (sky-blue curves); G9 (brown curves); a 
fc -5 / 3 power law is shown by the orange-dashed line; NF (SF) 
stands for normal-fluid (superfluid). 


we have obtained from E’, for p = 1 to 6, i £ (n, s); each 
row of this Table has two lines; the first and second lines 
contain, respectively, the values of £” and . Table HIO 
(Supplemental Material) shows that Cp = Cp> f° r P = 1 
to 6, for the runs G1-G9, because of the mutual-friction- 
induced locking of the normal fluid and superfluid veloc¬ 
ities in the inertial range. 

Figures [2] (a) and (b) show, respectively, plots of 
and versus the order p; in these plots the orange line 
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FIG. 2. Plots versus order p of the multiscaling exponents: (a) 
(p and (b) for the shell-model runs G3 (purple curve), G4 
(green curve), G5 (sky-blue curve), G6 (brown curve), and G7 
(magenta curve). PG1, PG2 are the classical-fluid-turbulence 
runs (yellow curve) and £ K41 = p/3 is denoted by the orange 
line. 


is the K41 prediction C* 41 = p /3 and the yellow line 
shows the multiscaling exponents Cp of classical (super¬ 
script c), 3D-fluid turbulence. The nrultiscaling expo¬ 
nents Cp, i £ (n, s), which we determine from the 3D- 
HVBK shell-model, show deviations from these devi¬ 
ations depend on the values of p n /p and B. Moreover, 
for the run G4 (p n /p = 0.4, B = 0.9838), the Q’s (green 
lines in Figs. [2] (a) and (b)) are close to C* 41 = p/3. For 
the run G3 (p n /p = 0.25, B = 1.08), the Q's (purple lines 
in Figs. [2] (a) and (b)) lie roughly between C* 41 and Cp! 
for the runs G5-G7, the differences between Q and (p 41 
and (p depend on p. 

To understand the dependence of the multiscaling ex¬ 
ponents Q, i £ ( n,s ), on p n /p (which includes the vari¬ 
ation of B with temperature), we plot, in Fig. [3j Cp, for 
p = 1 to 6, versus p n jp from our runs G1-G9. Figure [3] 



FIG. 3. Plots of (p, for p = 1 to 6, versus p„/p, from our shell- 
model runs G1-G9. For the purpose of reference, we show the 
value of a classical-fluid-turbulence exponent £/, for order p, 
by a horizontal, dashed line; different colors indicate different 
values of the order p. The black, dot-dashed lines indicates 
Cp 41 = P/3- 


shows that, depending on the values of p n /p, the behav¬ 
ior of the exponents can be classified roughly into five 
regions I-V (demarcated by grey, dashed, vertical lines 
on the plot). Region I (p n /p < 0.1): The values of Cp 
are close to the classical-fluid-turbulence exponents Q. 
Region II (0.1 < p n /p < 0.3): Cp > Cp, for p > 3 and, 
for p = 1,2, Cp — Cp- Region III (0.3 < pjp < 0.4): 
Cp ^ Cp* 41 . Region IV (0.4 < p n /p < 0.65): Cp show 
significant deviations from both Cp and C* 41 - Region V 
(p n /p > 0.65): Cp show a tendency to move towards CO- 
We now examine the dependence of the multiscaling 
exponents Cp, i € (n, s) on p n /p 1 while keeping the co¬ 
efficient of mutual friction B = 1.5 fixed, in runs Bl- 
B19. These runs allow us to classify the behavior of (/. 
i £ ( n,s ), as a function of p n /P , more clearly than the 
runs G1-G9. In the Supplemental Material, in Table m 
we list the values of Cp, i £ (n, s), which we extract 
from Sp (Eq. [20|, for p = 1 to 6, i £ (n, s); each row 
of this Table has two lines; the first and second lines 
contain the values of Cp and Cp, respectively. For these 
runs Cp — Cp- Fig- SI we plot Cp, versus p n /p, for 
p = 1 to 6 in runs B1-B19. These plots show two re¬ 
gions (0.1 < p n /P < 0.3 and 0.4 < p n /p < 0.65) with 
clear bumps, where the values of Cp deviate significantly 
from both C* 41 (< Cp ) and Cp(< Cp )• We classify roughly 
the behaviors of these Cp into six regions I-VI (demar¬ 
cated by grey, dashed, vertical lines in Fig. SI), which we 
describe below. Region I (p n /p < 0.1): Cp — Cp- R e_ 
gion II (0.1 < p n /p < 0.3): Cp differs significantly from 
both Cp and C* 41 , with Cp < Cp and C* 41 < Cp- Re¬ 
gion III (0.3 < p n /p < 0.4): Cp ^ C* 41 - Region IV 
(0.4 < pn/p < 0.65): Cp differs significantly from both 
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FIG. 4. Plots of (p, for p = 1 to 6, versus pn/p, from our 
shell-model runs B1-B19. For the purpose of reference, we 
show the value of a classical-fluid-turbulence exponent C^, for 
order p, by a horizontal, dashed line; different colors indicate 
different values of the order p. The black, dot-dashed lines 
indicates C^ 41 = p/3. In the shell-model runs B1-B19, we 
keep the mutual-friction coefficient B = 1.5 fixed. 


FIG. 5. Plot of (p, for p = 1 to 6, versus p„/p, from the shell- 
model runs R1-R12. For the purpose of reference, we show the 
value of a classical-fluid-turbulence exponent for order p, 
by a horizontal, dashed line; different colors indicate different 
values of the order p. The black, dot-dashed lines indicates 
Cp 41 = p/3. In the shell-model runs R1-R12, we keep the 
normal-fluid density fraction p n /p = 0.5 fixed. 


Cj, and C^ 4 \ with C p < Cp and C * 41 < <" Region V 
(0.65 > p n /p < 0.75): Q shows a tendency to move 
towards Region VI {p n /p > 0.75): Cp — Cp- 

We also explore the dependence of the multiscaling 
exponents Cp> * G (n, s), on the mutual-friction coeffi¬ 
cient B , while keeping the normal-fluid-density fraction 
Pn/p = 0.5 fixed. In our 3D-HVBK-shcll-model runs 
R1-R12, we systematically vary the values of B ; we list 
the values of Cp; * G (n, s) obtained from E^ (Eq. l20l) . 
for p = 1 to 6 , in Table El in the Supplemental Mate¬ 
rial; each row of this Table has two lines; the first and 
second lines contain, respectively, the values of Cp and 
Cp- In Fig. [5] we plot Cp versus B, for p = 1 to 6 , for 
the runs R1-R12; the exponents Cp deviate significantly 
from their classical-fluid-turbulence counterparts Cp; hi 
the range 1 < B < 3, with Cp > Cp; f° r P > 3, C" < Ci; 
and C 2 marginally larger than Cf- As B —» 0.1 (small 
values) and B —» 10 (large values) the multiscaling expo¬ 
nents Cp — Cp; because, in the limit B —» 0 , the normal 
fluid and superfluid are uncoupled; and for very large 
values of B, the coupling is so strong that single-fluid- 
turbulence results emerge. 

We have checked explicitly that all the values of Cp and 
Cp, which we have reported above, satisfy the convexity 
inequality Eq. m- We illustrate this in the plots of 
Fig. El 


IV. CONCLUSIONS 

We have carried out extensive numerical simulations 
of the 3D-HVBK shell-model, specifically to study the 




FIG. 6 . Plots of (a) f™ vx for the runs B1-B19 (B = 1.5); (b) 
//(,x for the runs R1-R12 (p n /p = 0.5), where // vx = ( P 3 — 
Pi)&P 2 ~ (P 3 -P 2 )& pi - (P 2 -Pi)C 2 P3 , i € (n, s), and we take 
pi = 1, P 2 = 2, and p 3 = 3. The multiscaling exponents Cp; 
i € (n, s), satisfy the convexity constraint, if // vx > 0 , for any 
three positive integers pi < P 2 < p 3 . The *-axis label in the 
above plots indicates the run index, e.g., Bl. 


multiscaling of structure functions in superfluid turbu¬ 
lence, because such multiscaling has been studied much 
less than its counterpart in classical-fluid turbulence. Ex¬ 
perimental investigations of turbulence in liquid helium, 
below the superfluid transition temperature T \, have pro¬ 
vided evidence for multiscaling, in the inertial range [ 6 IJ- 
[63]. These experiments have also motivated our study. 
Direct numerical simulations of models for superfluids, 
e.g., the Gross-Pitaevskii equation and the HVBK two- 
fluid equations, have not been able to cover the large 
range of length scales that are required to obtain reli¬ 
able data for high-order structure functions. Shell mod¬ 
els, based on the HVBK two-fluid equations, have been 
used to study the statistical properties of 3D superfluid 
turbulence in both 4 He [13,0] and 3 He-B [UHll; these 
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studies have elucidated the natures of energy spectra and 
fluxes, for both forced, statistically steady and decay¬ 
ing superfluid turbulence. The only detailed investiga¬ 
tion of the multiscaling behavior of structure functions is 
an HVBK-shell-model study [26} . This study has shown 
that, for p n /p < 0.1 and p n /p < 0.9, the multiscaling 
exponents are close to those in classical-fluid turbulence; 
whereas, in the range 0.25 < p n /p < 0.5, high-order 
mutliscaling exponents deviate significantly from, and 
are smaller than, their classical-fluid-turbulence counter¬ 
parts. 

Our extensive study of the 3D-HVBK shell model has 
shown that the nmltiscaling of structure functions in su- 
pcrfluid turbulence is more complex than that reported 
in R ef. J26l . Bowever, our results agree with those of 
Ref. [26| in that, for p n /p < 0.1 and p n /p > 0.75, 
the multiscaling exponents are close to the classical- 
fluid-turbulence values. Moreover, we find that there 
are two regions, with 0.1 < p n /p < 0.3 and 0.4 < 
p n /P < 0.65, where the multiscaling exponents are larger 
than their classical-fluid-turbulence and K41 counter¬ 
parts, i.e., > (p and Q > (^ 41 , i G (n,s). In the 


range 0.3 < p n /p < 0.4, these exponents are close to the 
K41 prediction, i.e., ~ Cp A1 - We have also investi¬ 

gated the dependence of the multiscaling exponents on 
the mutual-friction coefficient B , with p n /p = 0.5 fixed; 
our results show that, for small (weak-coupling limit) and 
large (strong-coupling limit) values of B, the multiscaling 
exponents tend to their classical-fluid-turbulence values, 
whereas, in the range 1 < B < 3, there are deviations 
from the classical-fluid-turbulence behavior > Q, for 
p > 3. We hope our extensive study of the multiscal¬ 
ing of structure functions in the 3D-BVBK shell-model 
will stimulate detailed experimental and DNS studies of 
such multiscaling in quantum-fluid turbulence in differ¬ 
ent quantum fluids. 
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Pn/p 

cr 

& 

G 

G 

G 

G 


(B) 

G 

G 

G 

G 

a 

Cl 

PG1 

— 

0.378 ± 0.004 

0.704 ± 0.006 

0.996 ± 0.009 

1.26 ±0.01 

1.51 ±0.02 

1.74 ±0.03 

PG2 

- 

0.383 ± 0.003 

0.714 ±0.005 

1.007 ±0.007 

1.27 ±0.01 

1.52 ±0.02 

1.75 ±0.02 

Gla 

0.0450 

0.378 ± 0.008 

0.70 ± 0.01 

0.99 ±0.02 

1.26 ± 0.03 

1.50 ±0.05 

1.73 ±0.07 


(1.5260) 

0.378 ± 0.008 

0.70 ± 0.01 

0.99 ±0.02 

1.26 ±0.03 

1.50 ±0.05 

1.73 ±0.07 

G1 

0.0450 

0.387 ± 0.003 

0.723 ± 0.004 

1.024 ±0.006 

1.300 ± 0.008 

1.56 ±0.01 

1.80 ±0.02 


(1.5260) 

0.384 ± 0.003 

0.721 ± 0.004 

1.022 ±0.006 

1.300 ± 0.008 

1.55 ±0.01 

1.80 ±0.02 

G2 

0.0998 

0.390 ± 0.003 

0.732 ± 0.005 

1.040 ±0.007 

1.32 ±0.01 

1.58 ±0.02 

1.83 ±0.03 


(1.3255) 

0.389 ±0.003 

0.731 ± 0.005 

1.040 ±0.007 

1.32 ±0.01 

1.58 ±0.02 

1.83 ±0.03 

G3 

0.2503 

0.372 ± 0.006 

0.71 ±0.01 

1.02 ±0.01 

1.32 ±0.02 

1.60 ± 0.02 

1.86 ±0.04 


(1.0765) 

0.372 ± 0.006 

0.71 ±0.01 

1.02 ±0.01 

1.32 ± 0.02 

1.60 ± 0.02 

1.86 ±0.04 

G4 

0.4004 

0.3309 ± 0.0001 

0.6663 ± 0.0001 

1.0046 ± 0.0001 

1.3446 ±0.0001 

1.6858 ± 0.0002 

2.0276 ± 0.0003 


(0.9838) 0.3310 ± 0.0001 

0.6664 ± 0.0001 

1.0044 ± 0.0001 

1.3441 ± 0.0001 

1.6847 ± 0.0002 

2.0259 ± 0.0002 

G5 

0.4994 

0.375 ± 0.004 

0.732 ± 0.005 

1.079 ±0.006 

1.418 ±0.008 

1.75 ±0.01 

2.06 ±0.02 


(0.9848) 

0.374 ± 0.003 

0.732 ± 0.005 

1.079 ±0.006 

1.417 ±0.008 

1.74 ±0.01 

2.06 ±0.02 

G6 

0.6003 

0.386 ± 0.003 

0.737 ± 0.005 

1.065 ±0.007 

1.37 ±0.01 

1.67 ±0.02 

1.94 ±0.03 


(1.0447) 

0.385 ± 0.003 

0.737 ± 0.005 

1.064 ±0.007 

1.37 ±0.01 

1.66 ± 0.02 

1.94 ±0.03 

G7 

0.6493 

0.384 ± 0.003 

0.734 ± 0.005 

1.060 ±0.006 

1.37 ±0.01 

1.66 ± 0.02 

1.94 ±0.03 


(1.1034) 

0.384 ± 0.003 

0.734 ± 0.004 

1.060 ±0.006 

1.37 ±0.01 

1.66 ± 0.02 

1.94 ±0.03 

G8 

0.6995 

0.383 ± 0.003 

0.728 ± 0.004 

1.047 ±0.007 

1.35 ±0.01 

1.63 ± 0.02 

1.90 ±0.03 


(1.1924) 

0.383 ± 0.003 

0.728 ± 0.004 

1.046 ±0.007 

1.35 ±0.01 

1.63 ± 0.02 

1.90 ±0.03 

G9 

0.7501 

0.381 ± 0.004 

0.718 ± 0.006 

1.027 ±0.008 

1.32 ±0.01 

1.59 ±0.02 

1.85 ±0.03 


(1.3267) 

0.380 ± 0.004 

0.718 ± 0.006 

1.027 ±0.008 

1.32 ±0.01 

1.59 ±0.02 

1.85 ±0.03 


TABLE III. Multiscaling exponents from our shell-model runs PG1, PG2, and G1 — G9; each row of the Table has two lines; 
the first and second lines contain, respectively, the values of ( p and C, p . In the second column, p n /p is the normal-fluid density 
fraction (first line) and B is the mutual-friction coefficient (second line, in parentheses). 
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Pn/P 


0 



C 2 



Cs 



Ci 



Cs 



Ce 


B1 

0.05 

0.387 

± 

0.002 

0.724 

"± 

0.004 

1.026 

± 

0.006 

1.303 

± 

0.009 

1.56 

± 

0.02 

1.80 

± 

0.03 



0.384 

± 

0.002 

0.720 

± 

0.004 

1.023 

± 

0.006 

1.301 

± 

0.009 

1.56 

± 

0.02 

1.80 

± 

0.03 

B2 

0.10 

0.389 

± 

0.003 

0.732 

± 

0.004 

1.042 

± 

0.006 

1.328 

± 

0.009 

1.59 

± 

0.02 

1.84 

± 

0.02 



0.388 

± 

0.003 

0.733 

± 

0.004 

1.045 

± 

0.005 

1.333 

± 

0.008 

1.60 

± 

0.01 

1.85 

± 

0.02 

B3 

0.15 

0.402 

± 

0.002 

0.766 

± 

0.004 

1.111 

± 

0.006 

1.44 

± 

0.01 

1.75 

± 

0.02 

2.05 

± 

0.03 



0.395 

± 

0.002 

0.761 

± 

0.004 

1.107 

± 

0.006 

1.44 

± 

0.01 

1.75 

± 

0.02 

2.05 

± 

0.03 

B4 

0.20 

0.380 

± 

0.003 

0.733 

± 

0.005 

1.071 

± 

0.006 

1.399 

± 

0.008 

1.72 

± 

0.01 

2.02 

± 

0.02 



0.376 

± 

0.003 

0.730 

± 

0.005 

1.069 

± 

0.006 

1.397 

± 

0.008 

1.71 

± 

0.01 

2.02 

± 

0.02 

B5 

0.025 

0.364 

± 

0.002 

0.711 

± 

0.004 

1.053 

± 

0.005 

1.392 

± 

0.009 

1.73 

± 

0.01 

2.05 

± 

0.02 



0.360 

± 

0.002 

0.708 

± 

0.004 

1.051 

± 

0.005 

1.390 

± 

0.007 

1.72 

± 

0.01 

2.05 

± 

0.02 

B6 

0.30 

0.334 

± 

0.002 

0.672 

± 

0.002 

1.011 

± 

0.002 

1.351 

± 

0.002 

1.690 

± 

0.003 

2.024 

:± 

: 0.09 



0.334 

± 

0.002 

0.671 

± 

0.002 

1.010 

± 

0.002 

1.350 

± 

0.002 

1.688 

± 

0.003 

2.01 : 

±1 

3.009 

B7 

0.35 

0.345 

± 

0.001 

0.683 

± 

0.002 

1.021 

± 

0.002 

1.359 

± 

0.003 

1.696 

± 

0.005 

2.03 

± 

0.01 



0.3385 

± 

0.0009 

0.677 

± 

0.001 

1.015 

± 

0.002 

1.355 

± 

0.003 

1.693 

± 

0.005 

2.02 

± 

0.01 

B8 

0.40 

0.340 

± 

0.002 

0.679 

± 

0.002 

1.019 

± 

0.003 

1.359 

± 

0.003 

1.695 

± 

0.006 

2.03 

± 

0.01 



0.339 

± 

0.001 

0.680 

± 

0.002 

1.021 

± 

0.002 

1.361 

± 

0.002 

1.699 

± 

0.006 

2.03 

± 

0.01 

B9 

0.45 

0.365 

± 

0.002 

0.712 

± 

0.002 

1.057 

± 

0.003 

1.397 

± 

0.006 

1.73 

± 

0.01 

2.05 

± 

0.02 



0.353 

± 

0.001 

0.699 

± 

0.002 

1.046 

± 

0.003 

1.389 

± 

0.005 

1.72 

± 

0.01 

2.05 

± 

0.03 

BIO 

0.50 

0.376 

± 

0.002 

0.727 

± 

0.002 

1.072 

± 

0.004 

1.410 

± 

0.007 

1.74 

± 

0.01 

2.06 

± 

0.03 



0.365 

± 

0.002 

0.720 

± 

0.002 

1.068 

± 

0.004 

1.408 

± 

0.007 

1.74 

± 

0.01 

2.05 

± 

0.03 

Bll 

0.55 

0.382 

± 

0.002 

0.734 

± 

0.003 

1.078 

± 

0.004 

1.415 

± 

0.006 

1.74 

± 

0.01 

2.05 

± 

0.03 



0.370 

± 

0.002 

0.727 

± 

0.003 

1.075 

± 

0.004 

1.414 

± 

0.007 

1.74 

± 

0.01 

2.06 

± 

0.03 

B12 

0.60 

0.380 

± 

0.003 

0.735 

± 

0.005 

1.074 

± 

0.008 

1.40 

± 

0.01 

1.72 

± 

0.02 

2.03 

± 

0.03 



0.379 

± 

0.003 

0.734 

± 

0.005 

1.073 

± 

0.008 

1.40 

± 

0.01 

1.72 

± 

0.02 

2.03 

± 

0.03 

B13 

0.65 

0.379 

± 

0.003 

0.724 

± 

0.004 

1.048 

± 

0.006 

1.36 

± 

0.01 

1.65 

± 

0.02 

1.93 

± 

0.03 



0.378 

± 

0.003 

0.723 

± 

0.004 

1.047 

± 

0.006 

1.36 

± 

0.01 

1.65 

± 

0.02 

1.93 

± 

0.03 

B14 

0.70 

0.375 

± 

0.004 

0.714 

± 

0.007 

1.02 

± 

0.01 

1.31 

± 

0.02 

1.58 

± 

0.03 

1.82 

± 

0.05 



0.375 

± 

0.005 

0.714 

± 

0.007 

1.02 

± 

0.01 

1.31 

± 

0.02 

1.58 

± 

0.03 

1.82 

± 

0.05 

B15 

0.75 

0.371 

± 

0.004 

0.706 

± 

0.006 

1.012 

± 

0.009 

1.30 

± 

0.02 

1.56 

± 

0.03 

1.80 

± 

0.04 



0.371 

± 

0.004 

0.705 

± 

0.006 

1.012 

± 

0.009 

1.30 

± 

0.02 

1.56 

± 

0.03 

1.80 

± 

0.04 

B16 

0.80 

0.370 

± 

0.004 

0.697 

± 

0.006 

0.998 

± 

0.009 

1.28 

± 

0.01 

1.54 

± 

0.02 

1.80 

± 

0.04 



0.369 

± 

0.004 

0.697 

± 

0.006 

0.998 

± 

0.009 

1.28 

± 

0.01 

1.54 

± 

0.02 

1.80 

± 

0.04 

B17 

0.85 

0.371 

± 

0.003 

0.698 

± 

0.005 

0.998 

± 

0.007 

1.28 

± 

0.01 

1.54 

± 

0.02 

1.79 

± 

0.03 



0.371 

± 

0.003 

0.698 

± 

0.005 

0.998 

± 

0.007 

1.28 

± 

0.01 

1.54 

± 

0.02 

1.79 

± 

0.03 

B18 

0.90 

0.372 

± 

0.004 

0.699 

± 

0.006 

0.996 

± 

0.008 

1.27 

± 

0.01 

1.53 

± 

0.02 

1.77 

± 

0.04 



0.372 

± 

0.004 

0.698 

± 

0.006 

0.996 

± 

0.008 

1.27 

± 

0.01 

1.53 

± 

0.02 

1.77 

± 

0.04 

B19 

0.95 

0.374 

± 

0.004 

0.700 

± 

0.006 

0.994 

± 

0.008 

1.266 

; =t 

: 0.01 

1.52 

± 

0.02 

1.75 

± 

0.03 



0.374 

± 

0.004 

0.699 

± 

0.006 

0.994 

± 

0.008 

1.265 

i ± 

: 0.01 

1.52 

± 

0.02 

1.75 

± 

0.03 


TABLE IV. Multiscaling exponents £ v from our shell-model runs B1-B19; each row of the Table has two lines; the first and 
second lines contain, respectively, the values of £p and £p. p n jp is the normal-fluid density fraction; we keep the mutual-friction 
coefficient B = 1.5 fixed. 
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B 


Ci 



Ci 



<3 



Ci 



Ci 



Ci 


R1 

0.10 

0.371 

± 

0.007 

0.71 

± 

0.01 

1.03 

± 

0.02 

1.33 

± 

0.03 

1.61 

± 

0.04 

1.89 

± 

0.06 



0.370 

± 

0.007 

0.71 

± 

0.01 

1.03 

± 

0.02 

1.33 

± 

0.03 

1.61 

± 

0.04 

1.89 

± 

0.06 

R2 

0.50 

0.366 

± 

0.007 

0.70S 

; ± 

0.01 

1.04 

± 

0.02 

1.36 

± 

0.02 

1.67 

± 

0.04 

1.98 

± 

0.05 



0.366 

± 

0.007 

0.708 

i=t 

0.01 

1.04 

± 

0.02 

1.36 

± 

0.02 

1.67 

± 

0.04 

1.98 

± 

0.05 

R3 

0.80 

0.348 

± 

0.002 

0.691 

± 

0.002 

1.034 

± 

0.003 

1.376 

± 

0.004 

1.715 

± 

0.007 

2.05 

± 

0.01 



0.347 

± 

0.002 

0.691 

± 

0.002 

1.033 

± 

0.003 

1.374 

± 

0.004 

1.712 

± 

0.007 

2.04 

± 

0.01 

R4 

1.00 

0.346 

± 

0.001 

0.688 

± 

0.002 

1.031 

± 

0.002 

1.370 

± 

0.004 

1.704 

± 

0.008 

2.03 

± 

0.02 



0.345 

± 

0.001 

0.688 

± 

0.002 

1.030 

± 

0.002 

1.368 

± 

0.004 

1.700 

± 

0.008 

2.02 

± 

0.02 

R5 

1.25 

0.361 

± 

0.002 

0.711 

± 

0.002 

1.057 

± 

0.004 

1.397 

± 

0.006 

1.73 

± 

0.01 

2.05 

± 

0.02 



0.359 

± 

0.002 

0.710 

± 

0.002 

1.056 

± 

0.004 

1.394 

± 

0.006 

1.72 

± 

0.01 

2.04 

± 

0.02 

R6 

1.50 

0.372 

± 

0.002 

0.727 

± 

0.003 

1.072 

± 

0.005 

1.408 

± 

0.008 

1.73 

± 

0.01 

2.04 

± 

0.03 



0.370 

± 

0.002 

0.725 

± 

0.003 

1.070 

± 

0.005 

1.405 

± 

0.008 

1.73 

± 

0.01 

2.04 

± 

0.03 

R7 

1.75 

0.384 

± 

0.002 

0.743 

± 

0.004 

1.086 

± 

0.006 

1.418 

± 

0.009 

1.74: 

±i 

3.002 

2.04 

± 

0.03 



0.382 

± 

0.002 

0.742 

± 

0.004 

1.086 

± 

0.006 

1.417 

± 

0.009 

1.73: 

±i 

3.002 

2.04 

± 

0.03 

R8 

2.00 

0.382 

± 

0.003 

0.738 

± 

0.004 

1.078 

± 

0.006 

1.404 

± 

0.009 

1.72 

± 

0.02 

2.01 

± 

0.03 



0.382 

± 

0.003 

0.738 

± 

0.004 

1.077 

± 

0.006 

1.403 

± 

0.009 

1.73 

± 

0.02 

2.01 

± 

0.03 

R9 

4.00 

0.376 

± 

0.004 

0.720 

± 

0.007 

1.04 

± 

0.01 

1.34 

± 

0.02 

1.61 

± 

0.03 

1.87 

± 

0.05 



0.376 

± 

0.004 

0.720 

± 

0.007 

1.04 

± 

0.01 

1.34 

± 

0.02 

1.61 

± 

0.03 

1.87 

± 

0.05 

RIO 

6.00 

0.368 

± 

0.004 

0.699 

± 

0.006 

1.008 

± 

0.008 

1.30 

± 

0.01 

1.58 

± 

0.02 

1.84 

± 

0.04 



0.368 

± 

0.004 

0.699 

± 

0.006 

1.008 

± 

0.008 

1.30 

± 

0.01 

1.58 

± 

0.02 

1.84 

± 

0.04 

Rll 

8.00 

0.370 

± 

0.004 

0.699 

± 

0.006 

1.000 

± 

0.009 

1.28 

± 

0.02 

1.54 

± 

0.02 

1.79 

± 

0.04 



0.370 

± 

0.004 

0.699 

± 

0.006 

1.000 

± 

0.009 

1.28 

± 

0.01 

1.54 

± 

0.02 

1.79 

± 

0.04 

R12 

10.0 

0.370 

± 

0.004 

0.696 

± 

0.006 

0.993 

± 

0.009 

1.27 

± 

0.01 

1.52 

± 

0.02 

1.76 

± 

0.04 



0.370 

± 

0.004 

0.696 

± 

0.006 

0.992 

± 

0.009 

1.27 

± 

0.01 

1.52 

± 

0.02 

1.76T0.04 


TABLE V. Multiscaling exponents ( p from our shell-model runs R1-R12; each row of the Table has two lines; the first and 
second lines contain, respectively, the values of ( p and ( p . B is the mutual-friction coefficient; we keep the normal-fluid-density 
fraction p„/p = 0.5 fixed. 



